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Abstract 



We have developed two methods to study an inverse 
scattering problem of electromagnetic waves in chiral media, 
for a spherical scatterer perfect conductor. The first 
called Far-field inverse problem was described analytical in 
our previous work. The second called Near-field method is 
given in detail in the present paper. It is a geometrical 
method based on the scattered fields. Using Near-field 
experiments, in which the scattered field is measured at the 
source, we solve the corresponding inverse scattering 
problem that is to determine the coordinates of the center 
and the radius of the spherical scatterer. 



Introduction 

We study an inverse scat- 
tering problem of electromag- 
netic waves with time harmo- 
nic dependence applying the 
method of Near-field [ 1 ] - [ 3 ] . 
Spherical electromagnetic wa- 
ves generated by a point 
source incident on a spheri- 
cal scatterer perfect condu- 
ctor in a chiral media. 

Knowing the incident and 
scattered wave fields in the 
inverse problem we are 
looking for the coordinates 
of the sphere center and its 
radius . 

The corresponding problem 
is solved by geometrical me- 
thods using either energy 
scattering cross-section, or 

FRIDAY 10 JANUARY 2014 Vl 



the scattered fields. In the 
first case we refer to the 
inverse Far-field problem [4] 
and in the second case to a 
Near-field inverse problem. 

The chiral materials exhi- 
bit the phenomenon of optical 
activity i.e., the phenomenon 
that in which the plane of 
polarizations of linearly po- 
larized light is rotated as 
the light passes through an 
optically active medium. 

The chirality is a proper- 
ty which is often found in 
nature and reflects the asym- 
metry in the spatial inver- 
sion. An electromagnetic wave 
into a chiral media is analy- 
zed in a counter-clockwise 

■19 FUNKTECHNIKPLUS # JOURNAL 



( LCP ) and a clockwise (RCP) 
Beltrami field. Using vector 
spherical harmonic functions 
and appropriate expansions of 
the Beltrami fields and by 
extension the development of 
spherical electromagnetic wa- 
ves in spherical wave fun- 
ctions, we calculate the 
exact solution of the scat- 
tering problem for a spheri- 
cal perfect conductor in chi- 
ral media, and the correspon- 
ding back scattering [3] -[4]. 

In the inverse problem we 
calculate the scattered Near 
-field, when the source is in 
position r = r 0 . A similar pro- 
blem for achiral materials 
has been studied in the works 
[ 1 ], [ 2 ]. 

Statement of the problem 

Consider a point source at 
the position r 0 that produces 
spherical electromagnetic wa- 
ves in a chiral media near a 
scatterer Q' , perfect con- 
ductor, i.e. the surface of 
the boundary condition is sa- 
tisfied 



A spherical incident ele- 
ctromagnetic wave (E“ c (r), 

H“ c (r)) with time harmonic de- 
pendence in accordance with 
the Bohren transformation, 
analyzed in spherical Beltra- 
mi fields Q*“(r) and Q^(r) 
as follows : 



Ef(r) = 0£(r) + <C e (r) 
CM = ^ (Qft(r) - Q£(r)) 



( 2 ) 



where q = (p/e) 1/2 is the in- 
trinsic impedance of the chi- 
ral medium. The Beltrami 
fields satisfy the equations, 
[5], [6], 



j Vx Q(.,r e (r) = Kl.Ql.,r 8 (r) 
[ Vx Q/?,r e (r) = -K R Qi., r 0 (r) 


(3) 


where y L , y R are wave 


numbers 


for Beltrami fields 
given by, 


and are 


k # k 


(4) 


Yi ~ l - fcp ' Yr ~ l + kfi 



n x (E inc (r) + E sc (r)) = 0 

( 1 ) 

r = a 

We consider that the exte- 
rior Q = R 3 /Q of the scat- 
terer is homogeneous chiral 
media with fixed chirality 
(3, dielectric constant e and 
magnetic permeability p . 



With /c = w(ep) 1/2 , co being the 
angular frequency. The indi- 
ces L and R denote the LCP 
and RCP fields respectively. 
The spherical incident Bel- 
trami fields with suitable 
normalization have the form, 
as defined in the following 
works: [4] (issue 2, p. 9, 

relations (4), (5)), and [3]. 
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If Qt":(r) and Q^(r) are 
Beltrami fields corresponding 
via transformation to Bohren 
in Ef nc and Ht nc , then the 

scattering problem for the 
perfect conductor can be 
formulated in the following 
way: be found Q™(r), Q™(r) , 
which belong to the space 
C 1 (Q + ) n C(Q + ) , such that: 



(E* c (r), H* c (r)) . The scattered 
electric field will be depen- 
ded on the polarizations p L , 
p R , (see [4] (issue 2, p. 
9, relation (6)) and will have 
the decomposition 

E r 0 C ( r I Pl' Pr) = Qm - 8 ( r I P L ' Pr) + 

sc - - ( 6 ) 

+ Q ,4 (r I P L , p R ) 




(iii) : 



f x Q*' r0 (r) + iQ (0 = 0(7) 
f x Q“ e (r) - iQ^(r) = o(i) 

r ->■ 00 



The limits on radiation 
conditions (5(iii)), are ta- 
ken uniformly in all dire- 
ctions r e S 2 , where S 2 is 

the unit sphere in M 3 and n 
is the outward normal unit 
vector perpendicular to the 
surface on the scatterer. 

The incident electromagne- 
tic wave (E^ c (r), H“ c (r)) on the 

scatterer Q generates the 
corresponding scattered field 



Where Q“ e (r | p L , p R ) and 

Q^, ( r I P L , P R ) are the corre- 
sponding scattered Beltrami 
fields which have the fol- 
lowing behavior, when r ^ 00, 
[5], [7], 

Q^r 0 ( r I Piv Pr) = K (V , A r )' 

... 1 ( ? ) 
•Skr.O'lPi, P R ) + 0(p) 

with A = L, R . The functions 
g L and g„ are the LCP and 
RCP far-field patterns respe- 
ctively, which are defined by 
the following relationship 
[7], 

g„,r.(f Ip,., p*) = 

jfix[y A V xE^(r'|p„p R )- (8) 

S 

- w A Y 2 E s r ^(r' | p L , pje r cfs(r') 
f-1, A = L 

With m A = \ , y 2 = y L y R . 

[ 1,A = R 
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Near-field Inverse problem 

Using the exact solution 
found in the paper [4] (issue 
2, pp. 10-12) and considering 
appropriate asymptotic forms 
of the Bessel and Hankel 
functions in low frequencies, 

i.e. |y„a|<scl, calculate asym- 
ptotic expressions of the 
coefficients a L n , a„ , b L n and 
b R n . Specifically we have [3], 



a L _ 1 + $k , .2 

n 2iZ 2 {2n + iy Yi ’ 



y L a 0 



R i (1 - Bk) n + 2 . x2 

fln ~~2< d + m n+A YLa) 



bL „ i (1 + P/c)- 
n 2nC n (1 - pk) n + 1 {Yr ’ 

y R a^Q (10) 

IjR i (- L_ P k) , n) 2n + 1 

K 2C(2n + 1) iVRa) 

Where 

<■„= l-3-5-(2 n - 1) = (2n)!/(2”n!) 

In the inverse problem we 
calculate the scattered Near- 
field in the source r = r 0 . By 
the relation (21), in work 
[4], with r = r 0 and for LCP 
incidence, we obtain [3], 



E r 0 c ( r a i Pl) = 

= Z n(n 2 + ~ B n a n{(^(wr a )x + h n {y L r B )ty) + i (h n (y L r B )x - b n (y L r 0 )ijj)} + (11) 

n= 1 

00 f~jf f] + 1) _ 

+ Z 2 B n a n{( h n(YR r e)x ~ h n {y R r B )\\>) + i (-fy,(yv e ) x - h n (y R r B ) ip) } 



Using the asymptotic rela- with y A r B ^>0, A = L, R and the 
tions: series. 



MM) Zn 

k 0 (y A r B ) ( y A rJ 



^ T 2n + 1 _ T_ 



k n (y A r B ) -n£ n 

MY A r B ) (y A r B y + 1 



( 12 ) 



2 T 2n + 1 r 3 (l + r 2 ) 



(1 - r 2 ) 3 






i(2n + l)r”*‘.iei^C) (14 , 

n-1 V- 1 - 1 ) 
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inpn + l (15) 



2/1 + 1 



Z 4-I I -I 

n 



t[t 2 In (1 - r 2 ) + 2r 2 - In (1 - r 2 )] 

“ r 2 - 1 

and because of the relation- 
ships : 



Yr 1 - P k 

V L 1 + ( 3 /c 



, 0 < r < 1 , 



where r = a/r 0 , we obtain from 
(11) [3], 

E r e C ( r 0 I PJ ~ 

(1 + m T 5 6 (17) 

2(/ L a) 2 (1 - r 2 ) 3 Pl 

Therefore the measure of the 
scattered field E* c (r 0 | pJ is 



E r 0 c (r G I P L ) 



, a t pn r . 

Pl 'l 2(v L a) II r f (18) 
y L a -+ 0 



Similarly for the RCP inci- 
dence, we have that 



Oo I Pr ) 



p)|, r 5 

PJ I 2(y R a) 2 (1 - r 2 ) 3 (19) 

y R a -> 0 



Finally, based on the rela- 
tionships (18) and (19), we 
conclude that 



p£.(r a I P j| ~ 

- (1 - li (20) 

2(y,a) 2 (1 - r 2 ) 3 1 ’ 

y A a — » 0, A = L, R 

Choose a Cartesian coordi- 
nate system Oxipz, and five 
point-source locations, name- 
ly 0(0, 0,0), A 1 (l, 0, 0) , A 2 (Q r 

1,0), and A,(0, 0, 21) , which 
are at (unknown) distances 
r 0 , r ± , r 2 , r 3 and r 4 , respe- 
ctively from the sphere's 
center K. The parameter 1 is 
a chosen fixed length. The 
sizes of the resulting five 

measurements |E^ r .(r J . | p A )| are 



M = — ^ E sc (rip) 

J (l-^p/c) 



Pi= r ~T' b = 



Where j = 0, 1, 2, 3, 4 . Therefo- 
re we have the following five 



measurements 




b 5 p i 
M = 3 

3 ( P 2 j-by 


(23) 


with Pj>b 2 >0. Applying the 


law of cosines in 
K0A 4 and using that 


triangle 



r ) = P) 1 



FRIDAY 10 JANUARY 2014 



Vl—23 



FUNKTECHNIKPLUS # JOURNAL 



N.M. BERKETIS 



we obtain 



P 4 = 2 + 2 P3-P0 ( 25 ) 

Because of (24), we take 

p e j - 3 pyb 4 + 3p 2 b 8 - 



a is too small so that a 1 , 
then the relationship (23) 
shows [3], 




so we have the following sy- 
stem of equations 



Thus in equation (26) iden- 
tified the p, with j = 0, 1, 2, 

3, 4 . If the radius a of the 
sphere is known, solve the where j = 0, 1, 2, 3, 4 . 
system of six algebraic equa- 
tions (25) and (26). Also if 
you consider that the radius 



p 4 = 2 + 2p 3 - p 0 
p 8 M - h 5 = 0 
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